The close-limit approximation to Neutron Star Collisions 



Gabrielle D. Allen (1) , Nils Andersson (2) , Kostas D. Kokkotas (3) , 
Pablo Laguna^, Jorge A. Pullin^ 4 ) and Johannes Ruoff' 5 ' 
t 1 ' Max Planck Institute for Gravitational Physics 
The Albert Einstein Institute, D-144^3 Potsdam, Germany 

( 2 ' Department of Mathematics 
University of Southampton, Southampton SOU IB J, UK 
^ Department of Physics 
Aristotle University of Thessaloniki, Thessaloniki 54006, Greece 
' 4 - ) Center for Gravitational Physics & Geometry 
Penn State University, University Park, PA 16802, USA 
' 5 - ) Institut fur Astronomie und Astrophysik 
Universitdt Tubingen, D-72076 Tubingen, Germany 

Q\ '. 

, We develop a close-limit approximation to the head-on collision of two neutron stars similar to 

1 that used to treat the merger of black hole binaries . This approximation can serve as a useful bench- 

t-H ' mark test for future fully nonlinear studies. For neutron star binaries, the close-limit approximation 

> | , involves assuming that the merged object can be approximated as a perturbed, stable neutron star 

during the ring-down phase of the coalescence. We introduce a prescription for the construction of 
^5 | initial data sets, discuss the physical plausibility of the various assumptions involved, and briefly 

investigate the character of the gravitational radiation produced during the merger. The numerical 
results show that several of the merged objects fluid pulsation modes are excited to a significant 
(N \ level. 
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I. INTRODUCTION 

m ■ 

o : 

^\ i The inspiral and subsequent merger of compact binaries is one of the most promising sources of gravitational radi- 
Q\ ■ ation for the new generation of interferometric detectors (LIGO, VIRGO, GEO600, TAMA). Once these instruments 
""q j come online, we hope to learn much about physics in strong gravitational fields and at super- nuclear densities. This 
q-i is an exciting prospect, but we need reliable theoretical models of the relevant gravitational- wave signals to dig them 
out of what is likely to be a significantly noisy data stream. Hence, a huge effort is presently being made to model 
(3JT), the late stages of binary inspiral, both for black holes and neutron stars. For the inspiral phase the post-Newtonian 
approximation scheme has been pushed to very high orders jl]], and the available results are now at a rather sat- 
isfactory level of accuracy. For the final merger of the two binary companions, our understanding is not yet at a 
comparable level. To completely model the collision of two compact objects, one must resort to numerical relativity 
and fully nonlinear simulations assuming no symmetries. This provides a great computational challenge, and truly 
reliable results may not be available for several years. 

In the last few years, our understanding of black hole collisions has improved considerably. An important reason 
for the recent advances is the parallel development of the numerical approach and approximate methods. In fact, the 
so-called close-limit approximation (in which the late stages of the merger is modeled by considering the two black 
holes as a single perturbed one) is in remarkable agreement with the numerical simulations We thus have a 

powerful benchmark test for the fully numerical schemes, which gives us some confidence in the physical picture that 
emerges. Of course, a considerable amount of work remains before the black-hole problem is studied in complete 
generality. Most importantly, the calculations must be generalized to include rotational effects || . 

Compared to the collision of two black holes, the merging of two neutron stars involves, in addition to geometrody- 
namics, a non-trivial amount of "dirty" physics associated with the stellar fluid and a strong magnetic field . In the 
initial stages of collision, a shock wave will be generated; the merged object will heat up dramatically; the electromag- 
netic fields of the two stars will become intertwined; and a strong burst of energy may result. Even though the details 
remain to be understood, the merger is believed to lead to (most of) the observed cosmological gamma-ray bursts [||. 
From a gravitational point of view, this is interesting since it means that gamma ray bursts should be accompanied 
by a burst of gravitational radiation. Thus, there is a possibility that future gravitational-wave observations will help 
shed light on the gamma-ray puzzle. Anyway, at the later stages of the merger, several avenues are open. If the 
merged object is too massive to form a neutron star, it will collapse to form a black hole. But for an interesting range 
of masses the collapse may be temporarily halted. A supermassive neutron star can be supported either by thermal 
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pressure or by rotation In these cases, the final collapse will occur only after the object either cools or spins 

down. The likely extreme rotation of the merged object is relevant also for other reasons. A rapidly rotating stellar 
object can be both dynamically and secularly unstable. In the first case, the star spins so fast that matter is ejected 
from the equator. Via the so-called bar- mode instability, the core will shed some of its mass and a sufficient amount 
of angular momentum for it to become dynamically stable H . Once sufficient angular momentum has been shed to 
bring the rotation rate of the merged object below the dynamical Kepler limit, secular instabilities may come into 
play Based on recent results, we would expect the instability associated with the merged objects r-modes to spin 
the core down to a period longer than 10 ms in a year or so [jl0-13|. The merged object should, in fact, evolve almost 
exactly like a newly born neutron star (up to the eventual collapse in the case of supermassive configurations). Hence, 



the formation of a massive neutron star via merger of two less massive ones could be followed by a detectable J_4 15 
gravitational wave signal due to the r-modes. 

Clearly, it is a very difficult task to model all the aspects of the merger in general relativity On the other hand, the 
payoff from detecting and understanding these events is enormous. Neutron star mergers provide unique laboratories 
for physics at an extreme level. Given an understanding of the single pieces of the merger puzzle, we can hope to put 
together an accurate picture of the entire event. This paper presents an investigation into the late stages of binary 
neutron star mergers. Inspired by the success of the close-limit approximation for colliding black holes, we develop 
an analogous framework for the case of neutron stars. Although the close-limit approach to compact object mergers 
is unlikely to be as successful in the case of neutron stars (see later discussion), our hope is to gain insight into the 
character of the gravitational waves that emerge from these events. Furthermore, the framework that we develop 
should serve as a useful benchmark test for future fully general relativistic hydrodynamical simulations. 

In the close-limit approach, the late stages of binary merger are studied by assuming that the initial configuration 
can be viewed as (in some sense small) perturbations of the final object. The background is either a known analytic or 
a simple numerical solution. For black hole binaries, the background is a Schwarzschild or Kerr spacetime, depending 
on whether the remnant is a slowly or a rapidly rotating black hole. For the situation considered in the present 
paper, head-on neutron star mergers, the background can be either a black hole or a non-rotating relativistic star. 
The approximation that we will develop is only relevant in the latter case. The case when a black hole is formed 
requires a study of gravitational collapse, and is beyond our present aims. One may argue that the collision of two 
neutron stars should typically lead to an object that is too massive to support itself against gravity, and which must 
therefore collapse to a black hole. If that were the case, we would be restricted to the presumably small subset of 
mergers that involve less massive neutron stars. However, as mentioned above, supermassive neutron stars can (at 
least temporarily) be supported by both thermal pressure and rotation. Since the merged object is likely to heat up 
to temperatures beyond 10 11 K and ought to spin close to the Kepler limit due to conservation of angular momentum, 
these effect will be highly relevant. If the object is too massive to form a cold non-spinning neutron star, collapse 
will inevitably follow either on the cooling timescale (of the order of 10 s for cooling due to neutrino emission) or the 
spin-down timescale (which depends entirely on the mechanism due to which the merged object spins down) . Since 
both of these time-scales are orders of magnitude longer than the dynamical timescale of the merged object (typically 
at the ms level) the immediate aftermath of merger ought to be attainable within our close-limit approach. 

In the case of colliding black holes, the recent studies show that the emerging gravitational waves are dominated 
by the quasinormal modes of the final black hole. Intuitively, we expect a similar result for merging neutron stars. 
The merged object will pulsate wildly. This then leads to a characteristic gravitational-wave signal containing several 
of the nonradial modes of oscillation of the newly formed star. Provided that these modes can be identified in the 
gravitational-wave data, their particulars can be used to put strong constraints on the mass and the radius (as well 
as the equation of state) of the star |l^-[l^| . This is an interesting idea, but it requires that the modes are excited to 
detectable levels. Previous simulations indicate only that the fluid /-modes and p-modes, as well as the gravitational- 
wave u>-modes are excited when a non- rotating star is perturbed generically [l9| , p0[ | . With the present approximation, 
we can estimate the actual level of excitation of these modes at the late stages of neutron star merger. This will 
provide information that can be used to assess the detectability of these events. 

Time evolutions of black hole and neutron star perturbations have recently received considerable attention. For 
non-rotating black holes, the dynamics of perturbations can be investigated via the Zerilli equation pl| . Studies of 
such problems have led to a much improved understanding of the dynamics of non-rotating black holes. A similar 
approach has also been used for perturbed rotating black holes. Rotating black-hole perturbations have been studied 
in the time-domain by solving the Teukolsky equation |p2| . For perturbations of non-rotating neutron stars, there are 
many possible formulations, even within a particular choice of gauge. The results we present for the evolution of initial 
data representing neutron star head-on mergers were obtained using two independent codes based on perturbative 
variables within the Regge- Wheeler gauge p!9j , p0| . 

Given a reliable numerical code to evolve perturbations of a single neutron star, the close-limit calculation only 
requires the construction of relevant initial data. An obvious fundamental difference between close-limit initial data 
sets for black hole collisions and those involving neutron stars is the presence of a horizon. In fact, the success of 
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the close-limit approach for black holes may to a large extent rest on the presence of the horizon. The reason is the 
following: The initial data for two colliding black holes corresponds to rather large perturbations only in the region 
inside the peak (roughly at r = 3M) of the curvature potential. In effect, the perturbative scheme is unlikely to be 
reliable in this region. However, the bulk of these perturbations never escape to infinity. They are scattered by the 
effective potential and subsequently swallowed by the final black hole. Thus, the anticipated errors in the scheme 
have a minute effect on the gravitational radiation that reaches the far zone. Unfortunately, the neutron star case 
is not masked in this way. Since one would typically only include dissipation due to gravitational radiation in the 
simplest models, all "errors" in our initial data must escape through this channel in the absence of a central horizon. 
Consequently, the outcome of close-limit neutron star collisions are likely to be rather sensitive to the particular 
structure of the initial data. It should also be pointed out that, the fundamental qualitative difference between the 
two problems makes a direct comparison between the black hole case and that of neutron stars rather dubious. 

However, a more careful study of the black hole results suggests that the close-limit approach for neutron stars may 
not be a lost cause after all. In close-limit black-hole collisions, most of the energy radiated, if not all, comes from 
quadrupole (I = 2) quasi-normal mode (QNM) ringing. The only aspect that distinguishes close-limit initial data 
from any other data is the fixed amplitude of the ringing. If we assume that a similar situation arises for neutron 
stars, we should focus our attention on the mode-ringing and ignore the initial burst of radiation. This is natural since 
one would expect the initial burst to be strongly dependent on the characteristics of the initial data. In the case of 
stars, this strategy is associated with an unfortunate sacrifice. The early parts of the gravitational-wave signal from 



the merger may contain the rapidly damped w-modes [23 13]. In the proposed prescription, the inferred amplitude of 



these modes may well be unreliable. Actually, one would expect this to be the case on physical grounds. The problem 
of constructing "astrophysical" initial data in numerical relativity essentially boils down to specifying the amount of 
gravitational waves in the spacetime. In a neutron star spacetime, it would be natural for the initial gravitational 
waves to escape via the w-modes. Hence, an uncertainty in the amplitude of these modes is expected, given our poor 
understanding of the relevant astrophysical data. Taking the suggested attitude, we do not expect to be able to assess 
the astrophysical relevance of the w-modes. Instead, we should view the neutron star close-limit initial data (and the 
associated results) as a constraint on the amplitude of the longer lived fluid pulsation modes. 

II. CLOSE-LIMIT INITIAL DATA 
A. Stellar Models 

A neutron star model in general relativity is obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) system 
of equations: 

47rr> (1) 
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~T = -{Pa + Po) -j- > (3) 



where the spacetime metric is given by 

ds 2 = -e 2v dt 2 + e 2A dr 2 + r 2 d9 2 + r 2 sin 2 9 dip 2 . (4) 

Above, po and po are the density and pressure, respectively, and m represents the mass inside radius r. It is related 
to the metric function A by 

e- 2 ^!-^. (5) 
r 

All quantities are functions of the radial coordinate r only. The construction of initial data for the close-limit 
approximation involves two steps: First, one needs to solve the TOV equations and obtain stellar models for both 
the colliding stars and the final configuration. The second step consists of a suitable superposition of colliding stars 
followed by a subtraction of the background star. This yields the perturbations which are then the focus of the 
evolution. With this in mind, we have denoted background quantities of the star formed by the merger by the 
subscript and perturbations of these quantities by the subscript 1. The TOV solutions of the colliding stars will be 
denoted by an asterisk (*). For simplicity, we will only consider polytropic equations of state p = K p r , where K 
and r are the adiabatic constant and index, respectively. The adiabatic index T is related to the polytropic index n 
by T = 1 + 1/n. In the specific example provided later we use T = 2 (n = 1). 
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B. Initial Perturbations and Constraints 



The construction of astrophysically relevant initial data in numerical relativity is an outstanding problem not 
only from the mathematical point of view. In the present context, we would like the data to represent the initial 
configuration for the ring-down phase after a collision. Without the non-linear evolution that precedes the ring- 
down stage, the specification of initial data becomes a non-trivial "guess." In the 3+1 formulation H] of Einstein's 
equations, initial data consist of the spatial metric g^, the extrinsic curvature Kij and the matter fields. These initial 
data contain too much freedom. Only four out of the twelve components in (gij, Kij) are fixed by the constraints. 
The remaining pieces, including the matter fields, are freely specifiable and single out the particular situation under 
consideration. Furthermore, it is usually not obvious which four components in (gij, Kij) should be obtained from 
solving the constraints. 

The standard procedure to separate the freely specifiable data from that fixed by constraints in numerical relativity 
is York's conformal approach [|25|. In addition, this method provides a natural framework for "superposing" solutions. 
This is a clear advantage given our present project. For instance, given data (gij, Kij) for individual black holes, it is 
possible to "add" these solutions and solve the constraints to obtain self-consistent, fully non-linear, initial data |P6fl . A 
similar procedure can be applied to neutron star binaries. As mentioned in the previous section, our problem involves 
constructing perturbative initial data. Once the background matter and spacetime is specified, the construction of 
perturbative initial data consists of finding solutions that satisfy the Hamiltonian and momentum constraints to the 
order required by the perturbative expansion. 

As in the non-linear case, close-limit initial data is completely characterized, for both black holes and neutron stars, 
by the freely specifiable data (data not fixed by the constraints) and boundary conditions. That is, both the free data 
and boundary conditions select from the class of solutions of the constraints those dealing with the superposition of 
colliding objects. In this paper, we focus the discussion on simplest possible case: the close-limit head-on collisions 
of neutron stars that are initially at rest. This means that the initial extrinsic curvature and matter current density 
vanish to all orders. Therefore, the momentum constraints are identically satisfied, and one is only required to solve 
the Hamiltonian constraint. More general classes of neutron star initial data are discussed in p7| . For vanishing 
extrinsic curvature, the Hamiltonian constraint reads: 

R=\6np, (6) 

where R is the Ricci scalar constructed from gij and p is the matter density. 

Following York's method, the spatial metric g^ and matter density p are conformally transformed according to: 

9ij = 4 9ij (?) 
p = r S p, (8) 
where hats denote conformal quantities. With the above transformation, Eq. (JsJ) takes the form 

8A0-i?0 = -167T/50- 3 , (9) 

where A = V 2 ;V J and Vi denotes covariant differentiation associated with the conformal metric gij. 
We will assume at this point that metric perturbations only enter via the conformal factor. That is, 

0-0o + 0i (10) 
9ij = 9% ■ (11) 

However, it is important to realize that this is not a physically motivated choice. The reason for choosing data that 
is spatially conformally flat (i.e. g\^ — 0) is by no means a necessary condition for our procedure to work. In fact, for 
neutron stars the choice gjj — could lead to a suppression in the excitation of w-modes [|l9| . Still, the calculations 
are simplified considerably by this assumption, and it allows a better comparison with the black hole binary case 
since the initial data for multiple black holes is usually constructed assuming conformal flatness ^8j. At one level, 
the availability of non-linear evolutions will remove this uncertainty since the free data can then be determined from 
the outcome of these evolutions. However, there will still be a corresponding choice to make in the specification of 
initial data for the nonlinear phase. 

The conformal transformations (Q) and (||), together with the perturbative expansions (|l(]) and ([n]) yield: 

9% = <t>U% (12) 
g% =4 0^ (13) 

Po = o " 8 /5o (14) 
Pi = ^ Pi -8 0^01 /5 . (15) 
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With the above perturbative expansions, the Hamiltonian constraint takes the form 



I Ac/>i - [R < 



>7TpO' 



= -16tt 



(16) 



In summary, given the background (<fio , , po), close-limit initial perturbations of head-on collisions of neutron stars 
initially at rest consists of only two quantities: the perturbation of the conformal factor, </>i, and the perturbation of 
the conformal matter density, p\. In the following, we view the density perturbation p\ as free data, with <j>\ obtained 
by solving the linearized Hamiltonian constraint Eq. (16). Thus, the density perturbation p\ fully characterizes the 
collision, and we now turn to its specification. 



C. A Recipe for Close-limit Superposition of Stellar Models 



The key contribution of the present work is to provide a recipe for obtaining the density perturbation p\ from a 
suitable superposition of isolated neutron stars. To represent a "true" close-limit approximation, this superposition 
procedure must be such that the perturbations vanish as the separation between the stars vanish. In order to superpose 
two neutron stars and solve Eq. (|l^) , it is convenient to perform a coordinate transformation that brings the 3-metric 
(0) into the isotropic (conformally flat) form: 



ds 2 = <po(df 2 + f 2 dd 2 + f 2 sin 2 ( 
This is accomplished by setting the conformal factor to 



dip 2 ). 



1/2 



and transforming the radial-coordinate according to 



df 
dr 



With the metric given by (fl7|), Eq. (16) can be trivially rewritten as a radial elliptic equation: 



f 2 flf 



dr 



1(1 



h 6 7T po O 



= -2tt 



b 3 pi 



where we have used a standard decomposition in spherical harmonics: 



p(f, 6, ip) = po (r) + ^2 Pi (?) Y im (0, <p) 

I 111 

4>{fA<P) = Mf)+^2<h (r) Ylrn (6, if) . 



(17) 
(18) 

(19) 

(20) 

(21) 
(22) 



Irn 



To proceed, we need to provide a superposition procedure from which p\ can be obtained. To do this, we make 
the following observation: All of the black-hole close-limit initial data sets considered so far have had not only the 
appropriate zero-separation limit but also the correct infinite-separation limit. That is, as the separation of the black 
holes increases, the initial data converges to that of two isolated holes. A similar behavior is in principle desirable 
for neutron stars. However, the situation for neutron stars is more complicated because there are no simple relations 
between mass, radius and density as for black holes. Hence, we are forced to use a somewhat ad hoc prescription. We 
neglect tidal deformations etcetera and use the following superposition of density profiles of isolated neutron stars: 



pin = p*^ - a + p^ + 1) - p^p - c ) P ^ + c ) 



1/2 



(23) 



Here p* represents the conformally transformed density profile of the individual colliding neutron stars displaced a 
distance in conformal space. For simplicity we are assuming identical colliding stars. It is straightforward to verify 
that the superposition of densities ( p3| ) satisfies both the zero-separation and infinite-separation limits. 
Assuming a small displacement (i.e. imposing the close-limit condition), we find that 



p*(f i ± e) = p*(r>) ± ? v, p*(f») + ~ e i k % v k p.(f) 



(24) 
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Thus 



p = p* + \ (f V, h? + \ C V, % p* 



(25) 



so the conformal density perturbation is given by 

px = p-Po 
= P* - po + 



(26) 



We choose coordinates such that the centers of the colliding stars lie on the z-axis, so the off-set vector £ l is given by: 



= i, j cos6>, -I sinO, 



Substitution of fl27h into (E6J) yields 



px{f, 9, ip) = p*(r) - p (f) + - £ 2 



cos 



df 



n,( r) J + cos 2 9 -^P*{f) + sin 2 9 \-j^P*(r) 



df 



Making use of y/Air y 00 = 1 and tt Y20 + 5 = cos2 ®, Eq. ( p8| ) can be rewritten as: 

pi(f, 9, ip) =^2p 1 (r)Yi m (9, tp) 



(27) 



(28) 



(29) 



I in 
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Z j - n 

+ 15^ 



df 



df 



1 d 



4 d 



f df 



r df 
Y 2O (0, <p). 



Y oo (0, p) 



(30) 



From the above results, we deduce that the conformal density perturbation p\ has two contributions. One is a 
monopole part (m = 0, 1 = 0), 



pi = V47T 



1 



d 2 

df 2 



<> -r, + -e\^P*+[±P*) -~p* 



4 d 



(31) 



Since a monopole perturbation does not lead to gravitational radiation, we will not consider this contribution further. 
The second part is a quadrupole (m = 0, I = 2) perturbation given by the last term in Eq. (|30|): 



df 2 



df 



1 d 

f df 1 



(32) 



This is the dominant radiative source which should be used in Eq. (|2C|) to construct initial data. Notice that the / = 2 
perturbation ( |32| ) has the correct limit pi — > when | — ► independent of the stellar model used for p* . On the other 
hand, for the monopole perturbation (|^) to vanish in the limit £ — * we must have p* — > po We will now discuss a 
specific model in which this is the case. 



D. Denning "Realistic" Colliding Star Models 

As we already discussed in the introduction, the final stellar object that is formed by merger will be rather different 
from the initial stars. It will certainly be hotter and most likely spinning more rapidly since the angular momentum 
of the inspiral orbit must be conserved. In the case we consider here, that of head-on collision of two non-rotating 
neutron stars, we obviously need not worry about rotational effects. But we still we need to estimate the changes 
in the equation of state as the stellar material heats up during the merger. To do this we must speculate what the 
outcome of the collision may be, and specify a relation between the final background spacetime including po and the 
density distribution p* of the individual initial stars. 
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In order to have the correct zero- and infinite-separation limits, we must consider the relation between the masses 
of the colliding stars and the mass of the background star. A simplifying and to some extent reasonable condition 
is to assume that the mass lost during the collision is not significant, so the total mass is approximately conserved. 
Specifically, this implies that the total mass computed, in the physical space, from the background density po and 
the total mass obtained from the superposition of densities /?* in ( p3| ) are roughly the same. As we shall see, in the 
close-limit approximation this condition implies that properties of the colliding stellar models (for example, the radius 
of the stars) depend on the separation. 

TOV solutions for polytropic equations of state are parameterized by the central density p c , adiabatic constant 
K and adiabatic index T. We assume that the collision does not modify the adiabatic index (and set T = 2 for all 
models). This assumption is consistent with available numerical results |29|| . Therefore, the models used to specify 
differ only in central density and adiabatic constant. We want these parameters to reflect changes in the equation of 
state as the temperature increases, in such a way that we retain the final, much hotter, star at zero separation. Our 
method for constructing polytropic equations of state for the initial stars (in terms of p c and K) in relation to that 
of the final star is inspired by the recent arguments of Shapiro || for head-on collision of neutron stars from rest at 
infinity. 

The starting point is to notice that TOV solutions exhibit the following scalings: 

M( Pc ) = M(p c ) K n ' 2 (33) 
R( Pc ) = R(p c ) K n ' 2 (34) 
p c = p c K~ n , (35) 

with M and R the total mass and radius of the star. Above, tildes denote dimensionless quantities. To use these 
relations in the close-limit approach, we recall that we want our scheme to be valid in two limits: It should lead to 
the expected results both in the limit of zero and infinite separation. To achieve this we assume that the mass and 
radius of the colliding stars (index *) are related to the background star (index 0) by: 

M* = Mo f^ip) (36) 
R* = Ro fc^) , (37) 



where r\ is a monotonic function of the separation £ in physical space and represents the "overlap" of the two stars. 
Quite naturally, 77 ranges from to 1 with r](£ — ► 0) = and 7y(£ — > 00) = 1. This construction ensures that 
our approximation satisfies the desired limits. However, we have as yet no information about the function rj for 
intermediate separations. Such information could be obtained from fully numerical studies of merging stars. In the 
present study we make a natural, albeit quite arbitrary, choice. We assume that in the close-limit regime the overlap 
function is linear in the separation of the colliding stars and use 

= 4-- (38) 
Substitution of (|3^) and ( |37| ) into (^) and (||J), respectively, now yields 



(40) 

Since the ratio M/R is a monotonic function of p c , Eq. (|4C|) yields 

-M R* , s 

-^ = -^ = 1. 41) 
M R Q 

Hence, from (|3^), ([H]) and ([35|), one obtains the following relationships between adiabatic constants and central 
densities: 



M R V 2 / V^* 
Eq. ( p^ ) implies 

M* _ M 
i?„ Rq 
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K * =K °^ L T i y n (42) 

p\ = P° c (^) 2 • (43) 

In the limit 77 — > 1 , these scalings reduce to the case of two polytropic stars colliding head-on from rest at infinity [||. 
Here, we will use ( [42| ) and to relate the polytropic equations of state for the initial and the final stars, also for 
intermediate separations. 

It is important to point out that, if the amount of mass lost during the collision is not negligible conservation of 
mass is not given by Mq = 2 M» . As is clear from Eq. ( |36| ) , conservation of mass would be given by Mq = 2 M* only 
in the case of infinite separation. For the situation we consider, conservation of mass is demonstrated by comparing 
the mass Mq computed from po to that from p in (|2^). It is the mass obtained from the superposition rule (^3|) that 
takes correctly into account the double counting in the overlap of the density profiles of the colliding stars. 



III. NUMERICAL RESULTS 



Our numerical procedure to construct initial data consists of first solving, using a fourth-order Runge-Kutta inte- 
grator, the TOV equations (|l])-(||) for the background and colliding stars, together with the coordinate transformation 
(|l9|). These solutions are then conformally transformed and used to compute p\. Finally, Eq. (20) is solved for the 
perturbation of the conformal factorial as a boundary value problem using a standard tridiagonal solver p0| . The 
boundary conditions for tfii in Eq. (|20| ) are: regularity at the axis, 0i|f=o — 0, and asymptotic flatness at infinity, 

-§f<fil \f=oa = 0. 

The main aim of the present paper was to provide a prescription for close-limit initial data in the case of neutron 
stars. But, even though we will not discuss a large sample of numerical evolutions here, it is clearly appropriate 
to illustrate typical results obtained for the proposed neutron star close-limit approximation. We have used two 
independent numerical codes jl9],^0) to evolve the relevant perturbation quantities from the initial data obtained in 
the previous sections. These evolution codes have been well tested, and we have verified that they lead to identical 
results in the present case. 

In Figures we show a typical initial data and gravitational waves resulting from the close-limit approximation. 
This particular case pertains to a final stellar configuration with p° = 2.69 x 10 15 g/cm 3 and Kq = 100 km 2 . For these 
parameters, the mass and radius of the background star are Mq = 1.24 Mq and Rq — 9.0 km, respectively. The initial 
colliding stars, which are displaced a distance 0.1i?o from the center of mass, follow from p* = 2.98 x 10 15 g/cm 3 and 
if* = 90.25 km 2 and have Ai* = 1.17 Mq and i?o* = 8.58 km. Figure shows the perturbed conformal factor </>i and 
density p\. The character of the emerging gravitational waves is exactly what one would expect: A sharp initial burst 
followed by slowly damped oscillations. The long-lived oscillations are associated with the various fluid pulsation 
modes of the final configuration. That this is the case is clear from Figure ^, where we show the Fourier transform of 
the waveform in Figure 2. In Figures I and I we can also see that it is rather difficult to draw any conclusions about 
the relevance of the short-lived u>-modes from our sample evolution. While the modes are clearly present in the early 
part of the signal if we graph the variable S as defined by Allen et al |L9| , they are absent in the corresponding Zerilli 
function Z which should be a reliable measure of the emerging gravitational waves. However, at the present time this 
is perhaps neither here nor there: As already mentioned, we expect the assumption of conformally flat initial data to 
suppress the w-modes in a crucial way. Hence, we should not attempt to draw any conclusions about their potential 
astrophysical relevance from the present results. 

We can, however, meaningfully discuss the longer lived fluid modes. The various peaks in the spectrum shown in 
Figure |^ correspond directly to the fluid pulsation modes of the final configuration, the lowest frequency mode being 
the fundamental /-mode and the next one being the first of the pressure p-modes. This is interesting further evidence 
that these modes will be clearly excited whenever a neutron star is dynamically perturbed, which is highly relevant 
considering the recently devised method for inferring stellar parameters from detected gravitational waves carrying 
the mode-signature |L7 18|. However, one must still prove that these modes carry sufficient energy to be observable 
by the new generation of gravitational wave detectors. To investigate this issue, one should study a larger sample 
of close-limit evolutions and perhaps also attempt a comparison to the black-hole case. Work along these lines is 
currently in progress, and we hope to report on it soon. 
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IV. FINAL COMMENTS 



We have developed a framework for studying merging neutron stars using perturbation theory. Specifically, we have 
devised suitable initial data for the very late stages of binary merger, when the two merged stars can be considered as 
a final configuration plus perturbations. This "close-limit" approximation is analogous to the one that has provided 
surprisingly accurate results in the case of colliding black holes. However, as we have discussed in some detail, it 
is not straightforward to devise a similar approximation in the case of neutron stars. Our chosen scheme respects 
some of the required physical constraints. It has the correct limits at infinite and zero separation of the two stars. 
Furthermore, we have tried to model the changes in the equation of state brought about by the merger in a simple, 
but seemingly appropriate way. Still, it must be remembered that this is just a first step and that one could 
potentially refine the close-limit idea considerably. 

Even though it is clear that the close-limit approach to neutron star collisions has severe limitations (it will certainly 
never completely replace fully nonlinear general relativistic hydrodynamics simulations) we think that it can prove to 
be of considerable use. On a technical level, it should be rather straightforward to use our initial data sets, combined 
with the perturbation evolutions, as benchmark tests for fully nonlinear evolutions in numerical relativity taking 
account of the detailed fluid dynamics. This would in fact be a very useful test for a nonlinear simulation since it 
would ascertain that the detailed dynamics associated with the stars various pulsation modes could be resolved. Also, 
one could clearly use the perturbation equations to evolve any neutron star process at the late stages, thus saving 
valuable computing time. Furthermore, it seems likely that we can learn some important physics from our results. 
Evolutions from all the close-limit data sets that we have so far constructed show that the fluid pulsation modes of 
the final star are excited to a significant level. As far as the potential excitation of the gravitational w-modes are 
concerned, our present understanding is far from satisfactory, but as we have pointed out, an investigation of this 
issue likely requires a true understanding of "astrophysical" initial data and a relaxation of the standard assumption 
of conformal flatness. 

In conclusion, it is worth pointing out that our close-limit framework can be extended to more general situations. 
In the present study we chose to restrict ourselves to the head-on collision of two stars that are initially at rest. 
These assumptions can conveniently be relaxed to allow the stars to have initial momentum. A generalization to the 
physically relevant case of slow rotation also seems possible. Work in these directions is in progress. 
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FIG. 1. Perturbations of the conformal factor cj>i and density pi for a close-limit collision of equal-mass stars with con- 
formal separation £ = 0.1-Ro- The TOV parameters for the background and colliding stars are: p° c = 2.69 x 10 15 g/cm 3 and 
Kq — 100 km 2 . For these parameters, the mass and radius of the background star are M = 1.24 Mq and Ro = 9.0 km, respec- 
tively. The initial colliding stars, which are displaced a distance 0.1i?o from the center of mass, follow from p* = 2.98 x 10 15 g/cm 3 
and K* = 90.25 km 2 and have M» = 1.17 M and R * = 8.58 km. 
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FIG. 2. Snapshots of various functions illustrating for the evolution of the initial data shown in Fig. 1. We show the function 
S defined as in [|L9|, the Zerilli function and the perturbed conformal factor T = 4r0i/</> u . 
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FIG. 3. The Fourier transform of the waveforms for the function S and the Zerilli function Z shown in Fig. 2. The two 
sharp peaks correspond to the lowest frequency fluid pulsation modes of the final configuration. The first peak belongs to the 
fundamental /-mode while the next one is the first of the pressure p-modes. It is notable that w-modes seem to be present in 
S but not in Z. This effect is further discussed in the main text. 
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